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^\ ABSTRACT 

Q . In order to examine the far-IR/radio correlation at high redshift we have studied the 

| H ' Spitzer 70 /im/far- infrared (far-IR) properties of sub-mJy radio sources from the 13^ 

C/5 I XMM-Newton/ Chandra Deep Field by redshift and galaxy type: active galactic nucleus 

(AGN) or star forming galaxy (SFG). We directly detect 70 /im counterparts (at > 3(t 
significance) for 22.5% (92/408) of the radio sources, while for the rest we perform 
. stacking analysis by redshift and galaxy type. For the sources detected at 70 fim we 

^ ' find that the median and scatter of the observed flux density ratio, qyo, are similar 

I to previous results in the literature, but with a slight decrease in qyo towards higher 

T-H , redshifts. Of the radio sources detected at 70 /im 8/92 were already classified as AGN, 

OO ■ but two of which maybe SFGs. For the stacked sources we obtain a detection for the 

I SFGs at every redshift bin which implies they have mean flux densities a factor ^ 5 

\^ • below the original 70 /im detection limit. For the stacked AGN we obtain a detection 

' only in our highest redshift bin (1 ^ z 5) where we may be sampling hot dust 

On . associated with the AGN at rest- frame 12 — 35 /im. The combined observed mean value 

■ of qjo for the SFGs (detected and non-detected at 70 /xm) decreases gradually with 

^ I redshift, consistent with tracks derived from empirical spectral energy distributions 

• i-H . (SEDs) of local SFGs. Upon closer inspection and when comparing with tracks of 

' appropriate luminosity, the values of qyo broadly agree at low redshift. However, at 

H I z ~ 1, the observed qro (for ULIRGs) is 2cr below the value seen for local ULIRGs 

. . . 1 tracks, implying a difference in the SED between local and z ~ 1 ULIRGs. At higher 

redshifts, the convergence of the tracks and the higher uncertainties in qyo do not 
allow us to determine if such a difference persists. 
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1 INTRODUCTION 

The tight correlation over many orders of magnitude be- 
tween the far-infrared (far-IR) and radio luminosity of star 
forming galaxies ( SFGs) has been w e ll studied in the lo- 
cal Universe (e .g. iHelou et al] Il985l : ICondon et al] Il99ll : 
lYun et al.|[200lh . The radio emission of normal galaxies is 
dominated by synchrotron radiation from r elativistic elec - 
trons and free- free emission from HIT regions (|Condonll 19921 ). 
Both these mechanisms are related to the presence of young. 
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massive (M > 8Mq) stars; radio non-thermal emission 
arises from electrons accelerated by supernovae from these 
stars and thermal emission from ionized HII regions. The 
mid- and far-IR emission arises from dust absorption and 
subsequent re-radiation of UV/optical li ght. Deviations from 
this correlation are s een within galaxies ( Hughes et al.ll200^ : 
iMurphv et al.l |2006| ). but these probably derive from the 
longer mean free path of the relativistic electrons com- 
pared to the dust heating UV photons and the scatter in 
the relation within local gala xies is similar t o that seen 
betwee n local galax i es (e . g. [Yuu et al.] l200ll) . Some au- 
thors jHughes et all l2006l : iMurphv et all l2006h have sug- 
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gested that the far-IR/radio correlation may also have some 
dependence on star formation rate (SFR) and SFR density 
within a galaxy. 

The far-IR/radio correlation is important a s it is used 
to define the radio luminosity/SFR relation fe.g. lBel3l2003l ) 
and can also allow the selection o f radio-loud active galactic 
nuclei (AGN) fe.g. iDonlev et all 2005). Hence, determining 
whether this relation deviates significantly or becomes less 
tight at higher redshifts and higher luminosities is of great 
consequence. In particular the interpretation of data from 
the next generation of radio surveys, which will find Milky 
Way-like star forming galaxies out to high redshifts, will 
rely on the application of the radio luminosity/SFR rela- 
tion. In the distant Universe the relation between the mid- 
IR and radio l uminosities has been shown to appr oximately 
hold to z ~ 1 l|Garrettll2003 : lAppleton et al.ll2004l ') based on 
ISO and Spitzer observations, although there are some dis- 
agreements about the proportionality constant . At z ^ 1 us- 
ing m illime tre and sub-mil l imetr e observations iKovacs et al] 
l|2006l ) and lVlahakis et al.l l|2007f) find evidence of a decrease 
of ~ 0.2 dex in the far-IR/radi o ratio compar e d to t he lo- 
cal value, but it was shown in ISevmour et al.l (|2008l . here- 
after SOS) that the relation does approximately hold by 
comparing AGN-free distant and local samples of powerful 
starbursts at the same rest-frame far-IR wavelengths. How- 
ever the SOS study was small consisting of eleven sources at 
1.5 ^ 2 ^ 3. 

Several different mechanisms could lead to departures 
from the radio/far-IR correlation at high redshift. D evia- 
tions at low luminosities have been noticed (jBelllbOOSi l . but 
these occur in low SFR sources (SFR^ 3MQyr~^) which 
cannot be detected to any great distance in current radio 
surveys. Modeling of the radio/far-IR relation in dusty star- 
bursts has shown that the correlation is a natural result 
in any starburst if synchrotron emission dominates inverse 
Compton, and the electron cooling time is sh orter than the 
analo gous fading time of the supernovae rate l|Bressan et al.l 
|2002| ). However, deviations may be possible in both the early 
phase of a starburst, when the radio thermal component 
dominates the non-thermal component, and in the post- 
starburst phase, when the bulk of the non-thermal com- 
ponent originates from less massive stars. Other potential 
mechanisms that could lead to deviations from the relation 
include: (1) evolution in metallicity and dust properties, 
with accompanying changes in the far-IR spectral energy 
distributions (SEDs); (2) evolution in magnetic field proper- 
ties; (3) quenching from the Cosmic Microwave Background; 
or (4) the effects of hot intraclu ster gas in dense environ- 
ments fe.g. lMiller fc Owenll200ll ). 

The most sensitive infrared probe of high redshift 
star formation with current facilities is the 24 /im band 
of the Multiband Imaging Photometer for Spitzer (MIPS) 
instrument (Rickc ct al. 200^ ) on board the Spitzer Space 
Telescope (jWerner et all 120041 ). The cross-correlation of 
radio and 24 /xm data from many deep surveys has 
led to several evaluations of a mi d -IR/radio relation 
fBovle et al.' '200?; 'Be swick et al] l2008l : |Garn & Alexander' 
2009 ; Ibar ot al. 200S). So me of these authors ( Boyle ct al. 
20071 : [ Beswick et alj |2008| ) observe a change in the mid- 
IR/radio ratio at faint radio and/or 24 /^m flux densities, 
but as these authors discuss there are different populations 
being analysed depending on the selection of the parent pop- 



ulation. Also, studies using 24 /im observations suffer from 
the shift of this band with redshift to shorter wavelengths 
where it no longer tracks the cool dust directly associated 
with star formation and into a regime where features due to 
polycyclic aromatic hydrocarbons (PAHs) may dominate. 
Furthermore, it is possible that different processes (star for- 
mation or AGN) dominate in each of the radio and 24 fim 
bands. 

We aim to overcome some of these difficulties by us- 
ing a sample of sub-mjy radio sources that have been 
categorised by the physical process that dominates their 
radio emission: star formation or AGN activity. In SOS 
we used several radio-related discriminators (radio mor- 
phology, radio spectral index, observed 24/xm/radio and 
/f— band/radio flux density ratios) to separate these two 
populations in a very deep 1.4 GHz survey and confirmed the 
expected resul t (from m odeling of the radio source counts — 
iRowan-Robinson et al.l ll991: Hopkins 2004) that the SFG 
population dominates at the faintest radio flux densities, 
5'i.4GHz ^ 0.1 mjy, but with a significant, ~ |, contribution 
from AGN. 

In this paper we examine the far-IR/radio correlation 
of the sub-mJy radio population by galaxy type and redshift 
by studying their 70 /im emission. Studying the emission at 
70 fim has several advantages over that at 24 nm. This far-IR 
band is more comparable to the infrared wavelengths used 
in the early studies of this correlation (i. e., using the IRAS 
60 and 100 /im bands: iHelou et allflOSSl ). The 70 ^im band 
also samples the SED of local (U)LIRGs closer to its peak 
(60 /im ^ Arcst ^ 120 /tm) than the 24 /tm band and is a 
more direct measure of the total IR luminosity and hence 
SFR. Furthermore, given the redshift range of the sample 
examined here, an analysis at 24 /xm would probe rest-frame 
wavelengths 7 /im ^ Arest ^ 24 /tm, well away from the cold 
dust related to the ongoing star formation, and into a regime 
containing uncertain tracers of star formation (e.g. PAHs) as 
well as complications due to silicate absorption/emission fea- 
tures and possible AGN contamination of unknown strength. 
The 70 /im band has a few disadvantages: its relative shal- 
lowness compared to equivalent surveys at 24 /im and the 
potential confusion at faint flux densities due to the low res- 
olution at 70 /tm. We can overcome the flrst of these issues 
through the use of stacking techniques on sources not de- 
tected at 70 /im. The latter problem is mostly overcome by 
the accurate (sub-arcsec) positions of our radio sources, the 
parent population we are investigating. 

We have chosen to adopt in this paper a similar philos- 
ophy to that used in SOS, where we make no assumptions 
about the SED of the sources we observe, but use observed 
flux density ratios as a function of redshift at all times, and 
then compare these to models. This approach is the most 
appropriate since at each redshift slice we are looking at 
sources with different luminosities which may not be directly 
comparable to each other. We present our data analysis and 
cross-correlation in §2 and in §3 we present the results of 
our cross-correlation and stacking analysis. We discuss our 
results in §4 and present our conclusions in §5. Through- 
out we use a concordance model of Universe expansion, 

= 1 - = 0. 3, ^0 = 1, and Ho = 70 kms"^ Mpc"^ 
(jSpergel et al.ll2003l ). 
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2 DATA ANALYSIS 

2.1 Radio Data and Separation by Galaxy Type 

Th e radio sample u sed in this work was originally presented 
in ISevmour et al.l (|2004l . hereafter S04) which described 
our deep Very Large Array (VLA) 1.4 GHz observations of 
the 1 3" XMM/Chandra Deep Survey Field (jMLHardy et al. 
1 19981 ). We obtained a rms noise of 7.5 ^Jy at the center 
of the single VLA pointing and found 449 sources at 4(t 
significance out to a diameter of SOarcmin. Over the last 
decade we have obtained deep follow-up optical, near-IR 
and Spitzer imaging. We have obtained spectra for 163/449 
radio sources and these spectra were used to calibrate our 
14-band optical/near-/mid-IR photometric redshifts for the 
radio sources (see Dwelly et al. in prep, for a summary). Our 
spectroscopy revealed one radio source to be a star which we 
remove from further consideration. Only sources that were 
detected in at least 4 bands (236 sources) have photomet- 
ric redshifts. Hence, our parent sample consists of 448 ra- 
dio sources of which 162 have spectroscopic redshifts, 236 
have photometric redshifts and 50 sources have unknown 
redshifts. 

In SOS we described the separation of our parent sam- 
ple into galaxy type: AGN and SFG. In that paper we made 
the assumption that one of these processes dominated the 
radio emission and that we could use diagnostics directly re- 
lated to the radio emission to separate these radio sources. 
These empirical diagnostics were radio morphology (from 
combined MERLIN/ VLA observations), radio spectral in- 
dex (from combined VLA 1.4/4.8 GHz data), and two ob- 
served flux density ratios plotted as function of redshift: 
24/im/radio and radio/T^— band. Our analysis in S08 found 
269/449 faint 1.4 GHz sources to be SFGs and the remainder 
to be AGN. 

Our radio catalogue is not complete to the 30 /iJy detec- 
tion limit due to radio instrumental effects. The attenuation 
of the primary beam of the VLA away from the pointing 
center is by far the strongest effect and causes a decrease 
in the sensitivity of the survey by a factor of ~ 2 at the 
edge of the field of view. In S04 we calculated correction 
factors due to this effect (as well as other more minor ones) 
and used these correction factors in determining the radio 
source counts; each source was given a weight corresponding 
to its detection plus a factor representing all the sources at 
the same fiux density, but not detected due to the effective 
decrease in the sky area at faint flux densities. These correc- 
tion factors were also used in determining the comoving SFR 
density history in SOS and we use these same factors again 
here. These factors only apply to sources below 100 piJy and 
are generally low, < 2 for 50 /iJy < Siaghz < 100 fiJy, but 
do reach values 2 at Siaghz < 50 /xJy. We discuss the 
impact of the use of these correction factors in section 13.51 



2.2 Spitzer MIPS 70 fim Data 

The MIPS observations of the 13" field consist of a deep scan 
map using all three bands (24, 70 and 160 ^m) which cov- 
ers approximately 0.5 deg x 1 deg. These observations were 
obtained in July 2005 as part of MIPS instrument team 
GTO time (PI G. Rieke, programme identification number 
SI). Our 1.4 GHz radio observations cover one end of this 



MIPS map. The 70 ^m Basic Cahbrated Data (BCDs) from 
the Spitzer Science Center were processed offline using the 
Germanium Reprocessing Tools (GERtQ), following the fil- 
tering techniqu es adopted for the extragalactic First Look 
burvey fxFLS: llYaver et aLlliooj ). In particular, negative 
sidelobes near bright sources were removed using a combina- 
tion of high pass time median and column filtering with the 
bright sources masked. The BCDs were then mosaiced in the 
standard fashion with MOPEX^ and the rms noise in the final 
70 /xm map was ~ 2 mjy per beam. A list of sources was pro- 
duced with the APEX^ software, which performs robust PSF 
fitting, down to 3cr (6mJy). Flux densities were divided by 
a factor of 0.918 to colour correct for a. v x Fi, = constant 
SED (the most appropriate for the sources we are detect- 
ing) as the MIPS absolute flux calibration references a 10'* K 
blackbody. The source list is incomplete at the 3 a level of 
significance and may not be free from spurious sources, but 
spurious sources are not a concern for our approach, where 
we simply intend to measure the 70 fj,m flux density of each 
member of our parent radio sample. We note that the proba- 
bility of an individual 70 /im source being matched by chance 
with a radio source is low (^ 1.5 %). 

Out of the 448 radio sources in the parent sample, 40 
were not covered by the 70 /xm image. From the remaining 
408 sources we find 92 with 70 /im counterparts within 5" 
from our 3 a catalogue. These radio sources with 70 fim 
counterparts all have redshifts and include eight AGN and 
84 SFGs. Hence, there remain 316 radio sources with only 
70 /im upper limits. 

2.3 Stacking of 70 ^m Data 

To push our analysis to fainter 70 fim flux densities, we made 
stacked postage stamp cutouts of the 70 /xm image at the lo- 
cations of radio sources that were not matched to sources 
in the 70 fim catalogue. Each cutout was 128 x 128 arcsec 
(32 X 32 pixels after resampling to 4 arcsec pixels in the 
mosaicing) in size. These cutouts were combined using a 
weighted mean (where the weight is the the inverse of the 
square of the local rms). 

The stacking was done in several bins of galaxy type 
(SFG, AGN) and redshift. In SOS we found redshifts for aU 
our SFGs so the sources with unknown redshifts comprise 
one bin at z = of 47 AGN which are too faint in the opti- 
cal and near-IR for photometric redshifts to be determined 
(none of which are detected at 70 ^m). The bins were ini- 
tially chosen to be of equal size in log(l-|-2), except some bins 
were amalgamated to ensure each bin had at least 30 sources: 
the two highest redshift SFG bins were combined and the 
flve AGN bins were combined into two. Such amalgamation 
is necessary in order to obtain a signiflcant improvement, i.e. 
^ \/30, in the sensitivity comp ared to the unstacked image. 
Similar to iHuvnh et al.l ()2007l ) , offset stacked images with 
the same number of sources per bin were generated by ran- 
domly choosing a nearby position (< 64" or < 3.5FWHM) 
in the 70 /im image for each stacked source. Two hundred 
randomly offset stacks were generated, and the uncertainty 
in the stacked flux density is taken to be the standard devi- 
ation of these 200 measured values. 

^ Available from |http : //ssc . spitzer . caltech.edu/postbcd/ 1 
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Figure 1. The observed far-IR/radio flux density ratio (q7o = log[57o/^m/>S'i.4GHz]) of our radio-selected population plotted as a function 
of redshift; sources with unknown redshifts are plotted at z = 0. Sources are separated by galaxy type (SFGs: red squares, AGN: blue 
open circles) and upper limits represent non-detections at 70 /im at a 3 cr limit of 6 mjy. Average individual uncertainties are indicated 
although we note the uncertainty in the redshift is considerably less for the ~ ^ sources with spectroscopic redshifts. The tracks of local 
(Ultra) Luminous InfraRed Galaxies, (U)LIRGs (derived from empirical observations, Rieke et al. 2009) are shown as reference. The 
ULIRG (log(LiR/L0) 12) has a higher SFR than the LIRG (log(LiH/LQ) ^ 11) by a factor of ten and has a higher observed value of 
qro until z ~ 2.5. We also show the track of a radio-quiet QSO from Elvis et al. (1994). 
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Figure 2. The observed qro as function of redshift for just the 
radio sources detected at 70 ^tm SFGs: red squares, AGN: blue 
open circles). The solid line indicates the best fit linear regression 
to these data (SFGs and AGN). We discuss the 70 ^tm detection 
of eight AGN in the text. 



3.1 Observed qro for the Whole Sample 

The IR/radio correlation (at mid- or far-IR wave lengths) is 
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commonly defined as qiR = log[5'iR/S'i.4GHz] (e^ 
200l| : lAppleton eral]|2004l : iBeswick etalllm 
2008). We examine our observed values of q7o in Fig. [T] as a 
function of redshift for all sources (including upper limits for 
the radio sources not detected at 70 yLxn). Fig.[T]figure shows 
a wide scatter for the whole parent population, although 
most values are upper limits due to non-detection at 70 /im. 
The sources with very low qyo upper limits, as seen for many 
of the AGN, are typically the brightest radio sources in our 
survey. As our radio survey is flux density limited, and hence 
suffers from some redshift/luminosity degeneracy, it is not 
appropriate to fit to qyo as a function of redshift without ac- 
counting for the different luminosity ranges probed at each 
redshift. We choose to plot tracks calculated from templates 
of local galaxies, but shifted to the corresponding redshift 
and convolved with the appropriate band-passes, as a di- 
rect comparison to our observed values of qyo. We overlay 
tracks derived from empirical local LIRG and ULIRG SEDs 
("Rieke et al.ll2009t) and an unobscured radio-quiet AGN SED 
(Elvis ct al. 1993)- The starburst templates are derived from 
the mean of many local galaxies of the appropriate luminos- 
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Figure 3. The observed value of qro plotted as a function of redshift and source type (SFGs: red squares, AGN: open blue circles) for 
70 fim detected so urces and stacked non-detections at 70 fj,m. Overlaid are star forming tracks covering a range of IR luminosities from 
iRieke et al. 1 1I2OO9I ') at luminosities indicated in the figure including those in Fig 1. These lines become thinner at the redshift a template 
of a given luminosity becomes undetectable in our radio survey. The thin grey lines represent tracks from templates at intermediate 
luminosities with interval of 0.25 dex. The mean flux density ratios of sources undetected at 70 /im, derived from the stacked 70 fim flux 
densities, are indicated by larger symbols with error bars (indicating the width of the redshift bin and the uncertainty in qro). The SFGs 
have a detected flux density in each redshift bin and hence a determinable value of q7o whereas AGN are only detected in the highest 
redshift bin. 



ity. We note an empirical M82 track would lie marginally 
below the LIRG track presented here. 



3.2 Observed qro for the Detected Sample 

The radio sources detected at 70 all have values of 
q70 > 1.5 (Fig. [2]) mainly due to the sensitivity limit of 
the 70 /xm observations. We calculate the median value of 
q70 for this subsample in order to compare to earlier work. 
We find a value of 070 = 2.13 ± 0.24 (using a biweight es- 
timator 



Beers et al.f|l990l) . very similar to that found by 



lAppleton et al.l (|200. 



q70 = 2.16 ±0.17 (covering a similar 
redshift range). Additionally we derive values separately for 
each galaxy type obtaining a similar value for the 84 SFGs 
(qro = 2.13 ± 0.24) and a slightly lower value for the eight 
AGN (q7o = 2.02 ± 0.22). While SFG and AGN resuhs are 
marginally statistically different the similarity in values is 
probably due to the depth of the 70 /xm data not being able 
to detect sources, regardless of type, with a low value of q7o. 
Naturally, the values of observed median q7o for SFG and 
AGN are only upper limits for the radio selected sample as 
a whole and the true scatter is certainly greater. We also fit 
these data against log(l -|- z) by a simple linear regression 
as this measure of distance more closely follows look-back 
time. We find for the whole sample and by galaxy type: 

q7o = 2.15 ±0.10- 0.72 ± 0.29 xlog(l + z) (all) (1) 

q7o = 2.14 ±0.10- 0.75 ± 0.32 xlog(l + z) (SFGs) (2) 

q7o = 2.31 ±0.21 - 1.59 ±0.52 X log(l±z) (AGN) (3) 



We observe a decrease in q7o with redshift for the whole 
of the 70 /im detected sample which is principally due to 
sources with very high values of q7o being found at low red- 
shift, but which become much rarer at higher redshifts. This 
fit is shown in Fig [21 The sample fitted here is naturally bi- 
ased against low values of q7o due to the 70 ^m detection 
limit. We find that the SFGs show a similar trend to the to- 
tal population, but although the AGN have a steeper slope 
they are almost consistent within the uncertainties and the 
slope is largely dependent on the lowest and highest redshift 
sources. All of these slopes are consistent with the decrease 
in observed q7o seen in local templates redshifted to earlier 
epochs, e.g. as seen in Figure O 



3.3 Observed q7() for the Stacked Sample 

From our stacked images we obtained flux densities ranging 
between 0.3 and 2.5 mjy with uncertainties of ~ 0.3 mjy. 
These flux densities correspond to 0.15 to 1.25 times the 
rms of the original image (2mJy) showing a typical factor 
of ~ 7 improvement in sensitivity. For each stacked measure 
of the 70 fim flux density we calculate a mean value of q7o 
by dividing by the mean radio flux density of the sources 
in a particular bin. In Fig. [3] we plot this value of q7o as a 
function of redshift for both individual sources and stacks. 
In the four SFG redshift bins we flnd signiflcant detections 
of 70 flux density which give values of q7o at the lower 
end of the distribution for SFGs. Hence, every SFG has a 
detection at 70 fj,m, either directly or indirectly. The stacked 
AGN are not detected in the lower (2 ^ 1) redshift bin nor 
in the unknown redshift bin, but they are detected in the 
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Figure 4. The mean value q7o for all SFGs plotted as a function of redshift where we have combined the ratio for individual sources 
with the stacked value for undetected sources weighting them by the number of sources (filled red squares). Open red squares indicate 
the mean value of q7o for our SFGs detected at 70 fj,m in the same redshift bins. Uncertainties are derived from measurement errors 
and Poisson statistics. The numbers beside each data point indicate the mean total IR, luminosity (in units of log[LQ] and summed 
over 8 — 1000 fim) derived from the radio SFRs and the Kennicutt (1998) relation. The black lines again represent the tracks from local 
templates from Rieke et al. (2009) at luminosities indicted in the figure. These lines become thinner at the redshift a given template 
becomes undetectable in our radio survey. The very thin grey lines represent tracks from templates at intermediate luminosities with 
intervals of 0.25 dex. The trend of a decrease in qyo toward higher redshifts is generally consistent with tracks using local SEDs. 



highest redshift bin (1 ^ 2: ^5 5). At each redshift however 
the stacked value of qro for AGN is at least 0.5 dex below 
that for SFGs. 



3.4 Combined Mean Flux Density Ratio 

In Fig. Owe note that the whole population of radio-selected 
SFGs including non-detections at 70 fim has a scatter, i.e. 
range of values, in qro of at least one order of magnitude 
below z = 1. There is also a global trend toward lower val- 
ues of qro at higher redshifts as seen with just the sources 
detected at 70 /im. The scatter for the AGN below z — 1 
is larger, at least two orders of magnitude. The track of 
qro for a local radio-loud QSO (|Elvis et al.lll993 ) remains 
very close to qro ~ —1 for this whole redshift range (off 
the plot for the figures presented here). We note that typ- 
ical values of q 7o for radio galaxi es are ^ —2 using data 
from both low |Dicken et al. 20081 ) and high redshift sam- 
ples (|Seymour et al.ll2007l )" 

Given that we can detect all the SFGs at 70 /xm, ei- 
ther individually or stacked, we can then examine the mean 
flux density ratio of the whole population as a function of 
redshift. We decrease the size of the redshift bin widths to 
half the prior value and make new stacked 70 /im images. 



but restrict our analysis to sources up to z = 3 as we have 
very few SFGs above this redshift. This decrease in bin size 
is important to provide fidelity in following any trend with 
redshift and is possible due to the larger number of sources 
per bin when including the detected and undetected sources. 
We combine the mean value of qyo from the detected sources 
together with the stacked value of q7o in these new bins by 
applying a weight scaled by the number of sources they rep- 
resent in each bin. We derive the mean redshifts of each bin 
with a similar weighting scheme in log(l -I- z) space. The 
results of this approach are shown in Fig. |4] which demon- 
strates a steady decrease in the observed far-IR/radio ratio 
toward higher redshift. This decrease is steeper than that 
seen for the de tected sourc e s. We overlay the tracks of lo- 
cal SFGs from iRieke et al.l (|2009l ) at different luminosities 
ilog[LiR/Le] = 10,11,12, 13). 



3.5 Selection Issues 

The work presented here, in particular Fig. |4j is naturally 
dependent on the initial separation of AGN and SFGs in 
SOS where a source was classified as an AGN by at least 
one of the methods discussed in the introduction. While we 
are confident that the selection is correct in at least a sta- 
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tistical sense it may fail on an individual basis as discussed 
in §4.1. Several of the discrimination methods in SOS (e.g. 
radio spectral index, radio morphology) were less effective 
for fainter sources, but this will change with up-coming data 
from surveys by e-MERLIN and LOFAR. Such data will be 
particularly important as SFGs are expected to dominate at 
fainter flux densities. 

One possible bias in our selection of AGN in this work is 
the use of the 24 /xm/radio flux density ratio as a diagnostic 
in SOS. We used this flux density ratio (q24) to select AGN 
as sources that were more than 5 a away from the lowest 
SFG track in q24-redshift space. This selection could poten- 
tially lead to a bias against SFGs with low values of q7o. 
We identify 55/179 AGN solely on their value of q24 (none 
of which were actually detected at 24 ^m and hence are up- 
per limits). When we repeated the whole analysis without 
the 24/im/radio diagnostic our mean values of qro for the 
SFGs in Figs El&ll decreased by 0.03 - 0.1 dex. As this ef- 
fect is small we conclude the q24 selection in SOS does not 
significantly affect the results presented here. 

As discussed in §2.1 our initial radio survey is not com- 
plete at the faintest flux densities due to the effective de- 
crease in survey area at the faintest levels. However, we have 
previously derived and used the necessary correction factors 
in S04. To asses whether their use here has any systematic 
effect we reran our code without these correction factors and 
found very little difference to our results. In Figs [3] & |4] we 
observed a decrease in q7o of ~ 0.05 dex for the highest red- 
shift bins with smaller differences at lower redshifts. Given 
the small size of this change and the successful use of these 
correction factors in previous publications we believe any 
uncertainties in these correction factors do not significantly 
affect our results. 



4 DISCUSSION 

4.1 Direct Detection of Radio-Selected AGN at 

70 ^im 

An interesting result in Fig.|3]is the detection of eight radio- 
selected AGN at 70 ^m. Given their redshifts, they have 
values of qro consistent with the SFG tracks. Hence, these 
sources are either AGN with unusually strong 70 ^im emis- 
sion or they are SFGs mis-identified as AGN in SOS. 

Looking more closely at the reasons for their original 
selection as SFGs in SOS we see that four were selected due 
to their very flat or steep radio spectral index, three due 
to their compact radio morphology and one due to both its 
low radio/24 ^m flux density ratio and its compact radio 
morphology. Of the sources selected by their spectral index, 
2/4 have values that are close to the accepted range for 
SFGs. Given the low significance of these particular sources 
(3 — 4 ct) at 4.8 GHz their spectral index is more uncertain, 
hence it is possible that these sources may in fact be SFGs. 
The other two sources fiagged by their spectral index have 
a very steep and a very fiat spectrum well away from the 
range of SFGs. The four sources selected due to their radio 
morphology have compact unresolved radio emission from a 
0.3" beam, but extended, elliptical optical morphology. The 
radio brightness temperatures of these sources are not high 
enough to definitely confirm an AGN origin to the radio 



emission and we cannot rule out the presence of compact 
nuclear starbursts. 

ft is also possible that some of these sources are hybrid 
AGN/SFG objects, i.e. they have a non-negligible contribu- 
tion from star formation to the radio (in SOS we made the 
assumption that one process wa s completely domina nt at 
radio wavelengths) or the IR fe.g. iMaiolino et al.lll995h . For 
example, a radio source with 30 % of its radio fiux density 
from star formation may still be selected as an AGN by some 
of the selection criteria in SOS, but the value of qro for the 
star formation component would be just 0.5 dex higher than 
the hybrid value, still within the scatter for SFGs. Hence, 
simple fiux density ratios, which currently are the dominant 
diagnostic at fainter fiux densities, are not capable of diag- 
nosing such hybrid objects. Given the short time-scale of the 
starburst and radio-loud AGN phase and the wide range of 
possible radio luminosities from both these processes, it is 
reasonable to expect the chances of detecting many hybrid 
sources in a blind radio survey to be low. 

We note that only 8/92 sources detected at 70 /im are 
potential hybrids and that two maybe SFGs. Given the 
statistical approach of SOS, to have a few potential mis- 
identifications is not unexpected. If we were to assign all 
these sources as SFGs our results in Fig. 3] would change 
by ^ 0.02 dex. The distribution of these 'hybrids' in the 
redshift /qyo plane is only marginally different from the other 
detected sources (they have slightly lower values of qro at 
a given redshift). However, only two have uncertain radio 
AGN indicators, and the rest likely could be genuine AGN 
or hybrids as there is some evidence that Seyferts and radio- 
quie t AGN do follow the far-IR/radio correlation (|Rov et al.l 
[1993). This complication indicates the importance of future 
high resolution and low frequency radio data in determining 
the relative contribution of SFG and AGN within a single 
object as well as for a fiux limited sample. 

4.2 Decrease of Observed qro at Higher Redshifts 

In Fig. U we see that at z = the qro values of the local 
empirical templates generally increase with total IR lumi- 
nosity in the 10^""^^ Lq luminosity range. This increase is 
mainly due to the more luminous starbursts having SEDs 
which peak at shorter wavelengths, i.e. they may be charac- 
terised as having hotter dust. Thi s trend is seen in empirical 
observations of l ocal starbursts (|Sanders fc Mirabe il Il996l : 
iRieke et al1l2009l ). The exception is in the 10^^"" Lq lumi- 
nosity range where some of the intermediate tracks in grey 
have higher values of qyo than the 10^^ Lq tracks while the 
qro for the 10^'^ Lq track is very close to the lO'^^ Lq track 
compared to these intermediate tracks. In this luminosity 
range local sources show a relatively fiat (i.e. a ^ —0.7, 
where Sv oa i>°') radio spectrum at frequencies below 7 GHz 
which is attributed to free-free absorption and is refiected in 
the most luminous galaxy templates used here. While our 
mean values of qro for the SFGs follow the general decrease 
of the SED tracks toward higher redshifts they mostly follow 
the LIRG (log[LiR/L0] ^ 11) track at low redshifts before 
moving to values in between the LIRG and 0.1 x LIRG track 
across 0.5 ^ 2 ^ 1.5. However this redshift range is where 
the inferred total IR luminosities from the radio or 70 fim 
fiux densities are in the high end LIRG and low end ULIRG 
{log[Li-R,/ Lq] ~ 12) regime. If our high redshift LIRGs and 



8 N. Seymour et al. 



ULIRGs had similar SEDs to their local counterparts we 
would expect values of qro higher by ~ 0.25 dex than we 
observe. 

This discrepancy between the data and the appropri- 
ate SED track by luminosity is a maximum of 2 cr in the 
z ~ 0.86 bin before the tracks begin to converge at higher 
redshifts. While we cannot completely rule out a few AGN 
contaminating our SFG sample we calculate that we would 
need at least 30% of the 52 SFGs in the z ~ 0.86 bin to be 
AGN with q7o ~ 1 to explain the observed difference. 

This difference in observed and expected values of qyo 
implies a change in the SED of z ~ f (U)LlRGs compared 
to local ones at moderate significance. For the ULIRGs one 
possibility is that there is less free-free absorption at high 
redshift which would increase their radio flux density rela- 
tive to the IR. This explanation could be true if star form- 
ing regions in high red shift ULIRGs are more extended (e.g. 
IChapman et al.l l2004f ) and hence less optically thick than 
local ULIRGs. By extrapolating the rest-frame 7 GHz un- 
absorbed radio power-law to the lower frequencies probed by 
our 1.4 GHz observations we can calculate the change in qro 
if no free-free absorption occurs in such objects. We flnd a 
decrease in the value of qro of 0.15(0.1) dex at z = 0.5(1.5). 
This explanation could account for a considerable amount 
of the difference in qro for the most luminous sources. 

Another possibility is that there is a change in the IR 
properties of LIRGs and ULIRGs at high redshift. There 
is evidence that a significant number of these sources at 
high redshift have different SEDs in comparison to local 
analogues, although the diffe rence is not well determined 
and depen ds on selecti on (e.g. Rowan- Robinson et al.|[2005l : 
IPapovicIi et al. 2007: Svmeonidis et al. 20081 ). For example. 
Symeonidis et al. (2008) present evidence that many sources 
selected at 70 ^im at 0.1 < z < 1.2 have SEDs that are pref- 
erentially cold, while Papovich et al. (2008) flnd that the 
SEDs of bright 24 /im sources at 1.5 < z < 2.5, on aver- 
age, are preferentially warm. It is likely that the structure 
of the luminous star forming regions are different in some 
(U)LIRGs at high redshift than locally, with more activity 
away from the compa ct nuclear sources tha t are typical for 
nearby ULIRGs (e.g. IChapman et al.l 120(3 ). This explana- 
tion would also permit the changes in radio emission dis- 
cussed in the previous paragraph. The work presented here 
is unable to distinguish between a change in the IR SED, 
a change in the radio spectral index or even some combi- 
nation of both. Determining the physical reasons for the 
differences between local and distant luminous starbursts 
will help us to understand the rapid change in star forma- 
tion rate density from the present day to 2 ~ 1 where these 
high luminosity objects d ominate star forming energy bud- 
get l|Le Floc'h et al.ir2005h . Future observations by Herschel 
will be able to characterise these differences in detail. 



4.3 Comparison with Other Studies 

Many of the rec ent studies of the mid and far-IR/radio 
correlation (e.g. iBovle et al ] I2OO7I : iBeswick et al ] l2008l : 
iGarn fc Alexander! 20091 ) have examined values of qiR as 
function of radio and /or IR flux density. Most IR/radio stud- 
ies have been based on 24 fim data because of its higher 
sensitivity, but as we discussed in the introduction we use 
70 /im data here as it more closely traces star formation. For 



our detected radio sou rces we found a media n value of qro 
comparable to that of I Appleton et all (|2004l ). but we find 
evidence for a decrease toward higher redshifts unlike the 
non k-corrected data from Appleton et al. (2004). We used 
data from a narrower, but deeper survey in the radio and 
mid-IR which could explain why we see this decrease. In the 
only other study to use 70 fim data Garn et al. (2009) exam- 
ined the radio fiuxes of 70 fim sources by stacking and found 
no evolution in qro over the 10 — 100 mjy 70 nm fiux range 
(similar to the fiux range of the 70 fim detected sources in 

this work) . 

llbar et al] l|2008l ) studied the 24/im/radio correlation 
by examining the 24 fim properties of a sub-mjy radio sam- 
ple. They found a decrease in observed (i.e. when no k- 
correction is applied) q24 with redshift, but not as steep 
as that seen here with qyo. The shallower slope is likely 
due to the SED of the selected sources being flatter at 
the wavelengths seen by the 24 /^m band. The 70 /im band 
used here probes longer wavelengths where the SED is ris- 
ing more steeply in typical starbursts, hence the observed 
value of q7o is likely to decrease more rapidly. Using an M82 
fc— correction those authors found an approximately con- 
stant value of q24 up to z ~ 3. This value, q24 = 0.71 ± 0.47, 
is lower than the value of qyo observed here at any redshift. 
This difference can be explained by the shorter wavelengths 
probed by the 24 /im band which are well away from the 
peak of the IR luminosity that traces the bulk of the star 
formation. The value of q24 also depends strongly on the 
SED assumed for the fc— correction as Ibar et al. discuss 
and, furthermore, this value of q24 is based on 24 pm detec- 
tions only and makes no other discrimination between AGN 
and SFGs. 

F or a sample of 24 /im selected sources, iBovle et al.l 
iiooj) observed a systematically higher value of q24 than 
seen by other authors. This value was constant with flux 
density. Given the selection wavelength and the lack of dis- 
crimination by galaxy type these authors are potentially 
selecting a large number of radio quiet AGN which have 
strong observed 24 /xm emission, but relatively little radio 
emissio n. However, in a sim ilar study of 24 selected 
sources iBeswick et al.l l|2008l ) found that the observed value 
of q24 decreased toward fainter 24 /xm flux densities and 
also decreased slightly toward higher redshif t s. In fact, at 
the faintest flux densities the IBeswick et al.l (l2008t ) result 
was almost an ord er of magnitude lower than that seen by 
iBovle et al.l l|2007l ). This discrepancy has not been fully ex- 
plained although it is possible the difference is a result of 
the different radio telescopes the different resolutions used 
(ATCA & MERLIN for Boyle et al. and Beswick et al. re- 
spectively). 



4.4 Comparison with Models of the Far-IR/Radio 
Relation 

While the ultimate origin of the far-IR/radio correlation is 
hot, young stars in dusty regions of the galaxy, the origin of 
the radiation we observe is not completely co-spatial or co- 
temporal. The wavelength range of far-IR radiation probed 
by original IRAS 60 and 100 fim band (and MIPS 70 fim 
band used here) originates from dust clouds relatively close 
to the young stars and emits fairly promptly after a localised 
burst of star formation whereas radio emission originates 
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from electrons traveling through the galaxy and interacting 
with its magnetic field which can happen relatively far and 
later than the original formation of stars. This difference in 
origin of the two wa velengths explain s the scatter we do see 
in the loc al relation jYun et al.ll200ll ) and the scatter within 
galaxies (|Murphv et al.ll2006l ). However, this difference also 
suggests that there is a dependence of the far-IR/radio cor- 
relation with the nature and mode of the very recent and 
inst antaneous star formation. 

iBressan et al.l (|2002l ) have performed detailed model- 
ing of the far-IR/radio correlation using the concept of age- 
selective obscuration (i.e. where younger stars and their cor- 
responding emission lie in denser molecular clouds) to derive 
SEDs and potential observables for a range of galaxy ages 
and SFRs. They find that the luminosity (far-IR or radio) to 
SFR ratio does vary by up to an order of magnitude with the 
age of the starburst for a range of star formation histories. 
Although these variations partly cancel out, these authors 
do predict a variation of qm, with starburst age which ap- 
pears stronger for starbursts with shorter e- folding times. 
They postulate that a low value of q/_R could be simply 
the natural consequence of a particular, e.g. post-starburst, 
phase in a galaxy's history. Hence, a further possible expla- 
nation for our low value of qro could s imply be the phase 
the starburst is in. iBressan et al.l l|2002l ) do predict an evo- 
lutionary path through the q/ii/radio spectral index param- 
eter space which could be used to determine the age of the 
starburst should sufficiently accurate IR and radio data be 
available in the future. 



5 CONCLUSIONS 

For the first time the far-IR/radio correlation of radio- 
selected SFGs and AGN has been examined to high red- 
shifts for a large number of sources. The 70 /im band is a 
more direct tracer of the bulk of the IR luminosity than 
24 /im and avoids many complicated spectral features in the 
mid-IR. Throughout this paper we adopted a philosophy of 
studying observed flux density ratios as a function of red- 
shift and comparing to the ratio observed from a redshifted 
SED rather than apply a fc— correction. This approach avoids 
the problem of choosing which template to use in making a 
fe— correction which can be problematic given the range of 
SEDs used previously for similar work, e.g. M82 and Arp220, 
and their application at high redshift. 
We have shown that: 

• Due to current radio and far-IR survey limits only 
sources with relatively high values of observed qro are de- 
tected at 70 fim. 

• The radio sources detected at 70 /im have similar ob- 
served median values and uncertainties in q7o as reported 
previ ously by a similarly selected sample (|Appleton et al.l 
I2OOJ), but also show evidence of a systematic decrease to- 
wards high redshift despite the bias towards high values of 
qro due to the 70 fim detection limit of the currently avail- 
able data. 

• Of those radio sources detected at 70 fim 84/92 were 
previously identified as radio-selected SFGs. On closer ex- 
amination due to their detection at 70 fim two of the remain- 
ing eight sources, identified as radio-selected AGN in SOS, 



are thought to be SFGs, or to show a significant contribu- 
tion from star formation, due to uncertainty in their radio 
spectra. 

• When stacking the sources not detected at 70 by 
galaxy type and redshift we find a stacked detection in each 
of the four redshift bins for the SFGs, but only detect the 
AGN in the highest redshift bin. In the lower redshift and 
unknown redshift AGN bins we do not find a stacked detec- 
tion down to ~ 0.3 mjy. The SFGs have stacked fiux densi- 
ties not much below our original 70 /im detection limit, hence 
the inferred values of qro are only slightly below those seen 
for the individually detected SFGs. 

• In order to examine the far-IR/radio correlation we de- 
termined the value of qro for all the radio selected SFGs. The 
mean values of qro for the detected and non-detected radio- 
selected SFGs were combined with appropriate weighting 
and we find an observed value of q7o that decreases toward 
higher redshift as expected from tracks derived from empir- 
ical SEDs of local SFGs. 

• When the observed values of qjo for the SFGs are com- 
pared closely to tracks of local SEDs at the appropriate lumi- 
nosity /SFR we find that while they broadly agree at low red- 
shift the observed qro at 2; ~ 1 (for sources close to ULIRG 
luminosities) is 2 a below the value seen for local ULIRGs 
shifted to this redshift. This result implies a difference in 
the SED of local and z ^ 1 ULIRGs, and demonstrates the 
value of future comparisons using bolometric luminosities 
when examining the far-IR/radio correlation at high red- 
shift. At redshifts higher than z ~ 1 the tracks of local SFGs 
converge and our observed values of qro have higher uncer- 
tainties, hence we are unable to determine whether such a 
difference persists to higher redshifts from this data set. 

We cannot completely rule out that part of this last 
result is due to AGN contaminants, although it would re- 
quire a large fraction of radio-excess sources. The lower 
value qro, beyond that simply inferred by the shifting of 
the band pass with redshift is consistent with other results 
in the literature that suggest that LIRGs and ULIRGs at 
high redshift are a more diverse population than their lo- 
cal counterparts (Rowan-Robinson et al."2005'; ' Saiina et al.l 
^2006; Brand ct al. 2008; Symconidis ot al. 200^. The fuU 
range of ULIRG properties at high redshift will be revealed 
by Herschel which will be the first telescope to truly select 
such objects based on their directly observed bolometric lu- 
minosities. 
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